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Abstract. Hirota bilinear method is proposed to directly construct periodic wave 
solutions in terms of Riemann theta functions for nonlinear equations. The asymptotic 
property of periodic waves are analyzed in detail. It is shown that well-known soliton 
solutions can be reduced from the periodic wave solutions. 

1. Introduction 

The investigation of the exact solutions of nonlinear equations plays an important 
role in the study of nonlinear physical phenomena. For example, the wave phenomena 
observed in fluid dynamics, plasma and elastic media are often modelled by the bell 
shaped sech solutions and the kink shaped tanh traveling wave solutions. The exact 
solution, if available, of those nonlinear equations facilitates the verification of numer- 
ical solvers and aids in the stability analysis of solutions. In the past decades, there 
has been significant progression in the development of these methods such as inverse 
scattering method [1,2], Darboux transformation [3-6], Hirota bilinear method [7-12], 
algebro-geometric method [13-18] and others. Among them, the algebro-geometric 
method is an analogue of inverse scattering transformation, which was first developed 
by Matveev, Its, Novikov, and Dubrovin et al. The method can derive an important 
class of exact solutions, which is called quasi-periodic or algebro-geometric solution, to 
many soliton equations such as KdV equation, sin-Gordon equation, and Schrodinger 
equation. In recent years, such solutions of nonlinear equations have been aroused 
much interest in the mathematical physics [19-24]. However, this method usually is 
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applied in the integrable nonlinear evolution equations admitting Lax pairs represen- 
tation and involves complicated algebraic geometry theory. These have been used far 
less than their soliton counterparts. The main reason for this is that they are far 
more complicated. Soliton solutions are typically expressed in terms of rational or 
hyperbolic functions, whereas qusi-periodic solutions require the use of Riemann theta 
functions and calculus on Riemann surfaces. Recently, Deconinck, Hoeij et al proposed 
an algorithm to compute the Riemann theta function and Riemann constants [25,26]. 

It is well known that the bilinear derivative method developed by Hirota is a 
powerful and direct approach to construct exact solution of nonlinear equations [7- 
9]. Once a nonlinear equation is written in bilinear forms by a dependent variable 
transformation, then multi-soliton solutions and rational solutions can be obtained. 
In recent years, Hirota method also has been developed for obtaining Wronskian and 
Pfaffian forms of N-soliton solution [10,11]. In this paper, we propose Hirota method 
for directly constructing periodic wave solutions in Riemann theta functions. It is 
shown that the periodic solutions can be reduced to classical soliton solutions under 
a certain limit. The appeal and success of this method lies in the fact we circumvent 
complicated algebro-geometric theory to directly get explicit periodic wave solutions. 
Moreover, all parameters appearing in the solutions are free variables, whereas usual 
quasi-periodic solutions involve Riemann constants which are difficult to be determined 
and need to make complicated Abel transformation on Riemann surface [13-18]. As 
illustrative examples, we then consider KdV equation and KP equation in this paper. 
Other many equations also can be dealt with this way. 

2. Periodic wave solutions of KdV equation and their reduction 

Consider KdV equation 



ut + Quu x + u xxx = 



0. 




Substituting transformation 



u = u + 2(\nf) 



XX 



(2.2) 



into Eq.(2.1) and integrating once again, we then get the following bilinear form 

G(D X , A)/ • / = (AA + Gu D 2 x + D x + c)f-f = (2.3) 

where c = c(t) is a integration constant and Mo is a constant solution of KdV equation. 
The Hirota bilinear differential operator is defined by 

D?D?f(x,t) ■ g(x,t) = (d x - d x ,) m (d t - d t ,) n f(x,t)g(x,t)\ x=x , tt=t , 

The D-operator have good property when acting on exponential functions 

D' x D r le^ ■ e 6 = (h - k 2 ) m {u; 1 - co 2 ) n e^ 2 

where £j = kjX + u>jt,j = 1, 2. More general, we have 

G(D X , A)e fl • e« 2 = G(k 1 - k 2 , u x - uj 2 )e^ 2 (2.4) 

We consider Riemann theta function solution of KdV equation 

j _ ^ ^ ^ni<rn,n>+2ni<!;,n> ^2 ^ 

ndZ N 

where n = (rii, • • • , n7v),£ = (£1, • • • ,£n), t is a symmetric matrix and Imr > 0, 
£j = kjx + ujjtj = I,-- - ,N. 

Consider the case when N — 1, then (2.5) becomes 

oo 



4 

Substituting (2.6) into (2.3) gives 

oo 

Gf ■ f = G(D x ,D t ) e 2 ™« + ™ v - e 



oo oo 

2 ^ ^ 2TTim^+Trim 2 t 



n=—oo m=— oo 

oo oo 



= E E G(D x ,D t y 

n=— oo m=— oo 

oo oo 

= E E G(27ri(n-m)A;,27ri(n-m)cj)e 27rl ( ri+m ^ +7r ' i (" 2+m2 ) T 

n=— oo m=— oo 

oo ( oo 

n+ "' \ E G(2«(2n-m0A;,2«(2n-m0cc;)e* + ^ m '^|>e 2 ^ 

m'=— oo ^n=— oo 

oo 

= ^ G(m')e Wf = 0. 



m — — oo 



Noticing that 

oo 

G(m')= G(2m(2n-m')k,2m(2n-m')uj)e^ n2+{n - m '^ 



n=n'+l 



n=— oo 
oo 



= +1 ^ G(27ri[2n'- (m' - 2)]k,2m[2n' - (m' - 2)]u) 

n'=— oc 

x exp(vri[(n /2 + (n' - (m' - 2)) 2 ]r) exp(2vri(m / - l)r) 

= G(m' - 2)e 2m(m '- 1)r 

G(0)e w(m '- 1) , m' is even 

( 5( 1 ) e ^(m'+i)( m '-2) ; m / ig odd 

which implies that if (5(0) = (5(1) = 0, then 

G{m') = 0, m f E Z 

In this way, we may let 

oo 

(5(0) = ^ (-167r 2 n 2 A;cj - 96M 7r 2 A; 2 n 2 + 256vr 4 n 4 A; 4 + c)e 2mr?T = 0, (2.7) 



n=— oo 

oo 



(5(1) = (-47r 2 (2n-l) 2 A;cu-24Mo7r 2 (2n-l) 2 A; 2 +167r 4 (2n-l) 4 A; 4 +c)e 7rl(2n2 - 2n+1)r = 0, 

?1= — oo 

(2.8) 



Denote 

S^n) = e 2 ™ 2r , 5 2 (n) = e M^-2n + i)r^ 

oo oo 

a n = — 16ir 2 n 2 k5i(n), a± 2 = 

n=— oo n=— oo 

oo oo 

61= (256vr 4 n 4 A; 4 -96Mo7r 2 A;V)(5i(n), a 21 = - ^ 4vr 2 (2n - lfA^n), 

n=— 00 n=— 00 

00 00 

a 22 = 6 2(n), b 2 = [!67r 4 (2n - 1) 4 A; 4 - 2Au ir 2 (2n - l) 2 k 2 }5 2 (n), 

n=— 00 n=— 00 

then (2.7) and (2.8) can be written as 

anu + a u c + bi = 0, 
a 2 i^ + a 2 2C + 6 2 = 0. 

Solving this system, we get 

6 2 ai2 - M22 M21 - b 2 a n 

uj = , c = . (2.9) 

^11022 — ai2a 2 i 011022 — o,i 2 a 2 i 

Finally we get periodic wave solution 

u = u + 2(lnf) xx , (2.10) 

where / and u are given by (2.6) and (2.9) respectively. From fig.l-fig.3, we see that 
the parameter r does affect the period and shape of wave, but parameter k has effect 
on the period and shape of wave. 

The well-known soliton solution of KdV equation can be obtained as limit of the 
periodic solution (2.10). For this purpose, we write / as 

/ = 1 + «( e 2 -« + e - 2 -«) + a 4 ( e 4 ^ + e" 4 ^) + • • • , 

where a = e tWT . 

Setting uq = 0, £ = £' — r/2, k! = 2mk, uj' = 27iiu, we get 

/ = 1 + e kWt + a 2 e~ 2 ^' + a 2 e 4m ^' + a 6 e~ 4 ^' + • • • 
— ► 1 + e k ' x+UJ '\ as a — >0 
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Figure 1. One-periodic wave for KdV equation and the effect of parameters on wave shape: 
(a) along x-axis, (b) along t-axis, where k = 0.1, r = i. 
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Figure 2. One-periodic wave for KdV equation and the effect of parameters on wave shape: 
(a) along x-axis, (b) along t-axis, where k — 0.1, r = 3i. 



So the periodic solution (2.10) convergence to well-known soliton solution 

U = 2(ln/) xx , / = l + e fe '* wt , 



only need to prove that 



a/ — > -k' 3 , as a — > 0. 



f2J 
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Figure 3. One-periodic wave for KdV equation and the effect of parameters on wave shape: 
(a) along x-axis, (b) along t-axis, where k — 0.15, t = i. 



In fact, it is easy to see that 

a u = -32k,7i 2 (a 2 + 4a 8 H ), 

a 12 = 1 + 2a 2 + • • • , 

a 2 2 = 2a + a 5 H , 

a 21 = -8A;7r 2 (a + 9a 5 H ), 

&! = 256A; 4 7r 4 (a 2 + 16a 8 + •••), 
6 2 = 16A; 4 7r 4 (2a + 81a 5 ), 

which lead to 

^2^12 — friO-22 = 32A; 4 7r 4 a + o(a), ana 2 2 — 012^21 = Skir 2 a + o(a) 
Therefore we have 

— > 4/c 3 7T 2 , as a — > 0, 



which implies (2.11). 



Now we consider two-periodic wave solution of KdV equation (N — 2). Substitut- 
ing (2.5) into (2.3), we have 

_ j _ ^ ^ (j(_D £)^g27ri<£,n>+7ri<T7i,n> . g27rj<5,m.>+7rj<rm,m> 

= G(27ri<n-m,A;>,27ri<n-m,cJ>)e 2m<f ' n+m>+m(<Tm ' m>+<T ™' n ^ 

oo 

n+ = =m ' ^ ^ G(27ri < 2n-m',k >,2ni < 2n-m',uj >) 

m'eZ 2 ni,n 2 =-oo 

x exp(7n(< r(n — m'),n — m! > + < rn, n >)) exp(27r? < £, m' >) 
= ^ GK.m^e 2 ^'" 1 ^ = 0. 

ro'GZ 2 

It is easy to calculate that 
G(mi,m' 2 ) 

oo 

= ^ G(27ri < 2n - m', k >, 2m < 2n - m', uj >) e ™«T(™-™'),n-m'>+<rn,n» 

ni,n2=— oo 

oo / 2 2 

£ G 2tu £(2nJ - (m; - 2<^) 2™ £(2nJ - {m' 3 - 2fy)) 

ni,n 2 =-oo \ j=l j=l 



UJj 



x exp I ni ^ [(nj + 5ji)T jk {n k + <J W ) + ((m^- - 2fy - raj) + S^T^m^ - 25 H - n' k ) + <J w j 

'(5(771; - 2, m ^ e 2^(m' 1 -l)r 11 +27rim 2 r 12) J = X 

G{m! x ,m! 2 - 2) e 27T ^ m 2- 1 ^+ 27Tim i T ^ , / = 2 
which implies that if 

(5(0, 0) = (5(0, 1) = (5(1, 0) = (5(1, 1) = 0, 
then (5(777^, m^) = and (2.2) is an exact solution of KdV equation. 



Denote 



dji = — 47r 2 < 2n — m? , k > (2n\ — m{)8j(n), 



lj2 



a j3 



i j4 



711,712=— oo 
oo 

An 2 <2n- m j , k > (2n 2 - m 3 2 )5j(n) 

711,712 = — OO 
OO 

24tt 2 < 2n - m j , k > 2 6j (n) , 

711,712 = — OO 
OO 

E 5 » 



711,712 = — OO 
OO 



bj = - ^ 4 <2n-m j ,k> i 5j{n), 

711,712 = — OO 

5 (fl) = f> 7Ti<T ( n ~ m ' i )' n ^ rn j >+iri<Tn,n> 

j = 1,2,3,4 m 1 = (0,0), m 2 = (l,0), m 3 = (0,l), m 4 = (1, 1) 
A = (a fci ) 4x4 , 6 = (6i, 6 2 , 63, & 4) T , 
then we have 



w 



b, 



from which we obtain 



Ai A 2 

LOl = — , W 2 



«0 



(2-12) 



A3 

A ' ~ z ~ A' " u ~ A ' 
where A = |A| and Ai,A 2 ,A 3 are produced from A by replacing its 1th, 2th, 3th 

column with b respectively. Finally we get two-periodic wave solution (See fig. 2) 



u = u + 2(\nf) xx , 



(2-13) 



where / and cji,cj 2 are given by (2.5) and (2.12) respectively. The properties of the 
solution are shown in fig.4-fig.7. Fig. 4 and fig. 5 shown effect of parameters k±,k 2 on 
period and shape of wave. Fig. 7 gives effect of parameter rn,r 12 ,r 22 on period and 
shape of wave. The fig. 7 and fig. 8 show that two-periodic wave can degenerate to 
one-periodic wave when k\ is sufficient small. 
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Figure 4. Two-periodic wave for KdV equation and the effect of parameters on wave shape: 
(a) along x-axis, (b) along t-axis, where th = O.li, t 12 = 0.2i, t 22 = 3i, k\ = 0.2, fc 2 = —0.3. 



(a) 



(6) 





Figure 5. Two-periodic wave for KdV equation and the effect of parameters on wave shape: 
(a) along x-axis, (b) along t-axis, where T\\ — O.li, t 12 = 0.2i, t 22 = 3i, fci = 0.1, ki = —0.3. 

The two-soliton solution of KdV equation can be obtained as limit of the periodic 
solution (2.10). We write / as 



_|_^ e 27ri(^i+^2) _|_ g-27ri(^i+^2)^ e 7r«(T-il+2ri2+r22) _)_... 



Setting = 2m^i — niT U , £' 2 = 2ni^ 2 — 7T2t 2 2, t 12 = ir (r is a real), we get 
— ► 1 + + e ? 2 + e ^+«2-2^ ) as aij a2 — > q 
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Figure 7. Two-periodic wave for KdV equation and the effect of parameters on wave shape: 
(a) along x-axis, (b) along t-axis, where m = 0.01, Tyi = 0.2i, r 2 2 = 3i, k\ = 0.1, k 2 = —0.2. 

where 

u'i—t—k®, as ai,ac2— > 0. 



3. Periodic wave solutions of KP equation and their reduction 

Consider KP equation 

U t = (U XX + 3U 2 ) X + W^Uyy. 



(3.1) 




Figure 8. Two-periodic wave for KdV equation and the effect of parameters on wave shape: 
(a) along x-axis, (b) along t-axis, where m = 0.001, T12 = 0.2i, t 2 2 = 3i, ki = 0.1, k 2 = —0.2. 



Substituting transformation 

u = u + 2(\nf) xx (3.2) 
into (3.1) and integrating once again, we then get the following bilinear form 

G(D X , D y , D t )f ■ f = (D x D t - 6u D 2 x - D x - 3D 2 y + c )f ■ f = (3.3) 

where c is a constant of integration and u is a constant solution of KP equation. Let 
us consider Riemann theta function solution of KP equation 

j ^ y ^ni<rn,n>+2ni<t;,n> ^ 

n€Z N 

where n = (ni, • • • ,Un),C, = (£1, • • • ,£n), t is a symmetric matrix and Imr > 0, 
= kj(x + p0 + Ujt) + = 1, • • • ,N. 

In the following, we consider some special cases. When N — 1, (3.4) is reduced to 

00 
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Substituting (3.5) into (3.3) gives 

oo oo 

Gf ■ f = G(D X , D v , D t ) e 27rm « +Trm2r • ^ ^mC+Wr 



OO OO 



= ^ ^ G(D x ,D y , £> t ) e 2™«+™« 2 T . e 2^ m ^m 2 T 

TL — — oo m— — oo 
oo oo 

= Z G(2^(n-m)fc,27r«(n-m)A:p,2^(n-m)Me 27rl(n+ " l)4+ " ( " 2+m2)T 

TL — — oo m — — oo 

oo f oo ^ 

n+,n=,n' ^ I G ( 2 ™(2n - m')fc, 2™(2n - m')fcp, 27ri(2n - m')Me™ [( " 2+( "~ m ' )2]r U 2 " m '« 

m ; — — oo ^n= — oo ) 

00 

= G(m')e 2W « = 0. 

m' — — oo 

Noticing that 

00 

G(m') = ^ G'(27ri(2n-mOA;,27ri(2n-m / )A;p,27ri(2n-m , )A;cj)e^ [(n2+{n - m ' )2]T 

n=— oo 

00 

n= = +1 ^ G(27ri[2n' - (m' - 2)]ib, 27ri[2n' - (m' - 2)]ifep, 2ni[2n' - (rri - 2)]ku) 

n'=— oo 

x exp(7ri[(n' 2 + (n' - (m' - 2)) 2 ]r) exp(27ri(m' - l)r) 

= G(m' - 2)e 2 " (m '- 1)r 

G(0)e w ( m ' _1 \ m' is even 

Q^ e ni(m'+l)(m'-2)^ m > ig odd 

which implies that if G(0) = G(l) = 0, then 

G(m') = 0, m' G Z 

In this way, we may let 

00 

(5(0) = ^ {-167r 2 n 2 k 2 u + 96u n 2 k 2 n 2 - 256n 4 n 4 k 4 + A87t 2 k 2 p 2 n 2 + c )e 2 ™ V = 0, 

n=— oo 

(3.6) 

00 

0(1) = (-4vr 2 (2n - l) 2 k 2 u + 2Au 7T 2 {2n - l) 2 k 2 - 167r 4 (2n - 1) V 

n=— 00 

+ 12vr 2 (2n - l) 2 A; 2 p 2 + Co ) e ^(2« 2 -2n+i)r = Q? (3^ 
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Denote 



oo 

2n 2 



a = e nrt , a u = -16tt 2 A; 2 n 2 a 2n \ a 12 = J] 

n=—oo n=—o 

oo 

& i = - E (96^ 2 n 2 fc 2 - 256 7 rVA : 4 +48^ 2 fcV« 2 )« 2 " 2 , «2i - -47r 2 fc 2 £ (2n-l) 2 a 



2n -2n+l 



a 22 = Yl a2n2 ~ 2n+1 > b 2 = E [167r 4 (2n-l) 4 fc 4 -24 U0 7r 2 (2n-l) 2 /fc 2 -12^ 2 (2n-l) 2 fc 2 p 2 ] a 2 " 2 - 2n+1 , 

n— — oo n— — oo 

then (3.6) and (3.7) can be written as 

a u uj + a 12 c = &i, 
a 2 iu; + a 2 2C = &2- 

Solving this system, we get 

M22 - b 2 a 12 b 2 a n - b x a 2X 
lo = , c = . (3.8) 

a ll a 22 — a 12 a 21 a ll a 22 — a 12 a 21 

Finally we get periodic wave solution of KP equation (3.1) 

u = u + 2(lnf) xx , (3.9) 

where / and u are given by (3.4) and (3.8) respectively. From fig.9-fig.ll, we see that 
the parameter r does affect the period and shape of wave, but parameter k has effect 
on the period and shape of wave. 

The well-known soliton solution of the KP equation can be obtained as limit of the 
periodic solution (3.9). To this end, we write / as 

/ = 1+ a (e 2m t + e" 2 ^) + a 4 (e 4 ^ + e~ 4 ™«) + • • • . 

Setting uq — 0, 7 = 70 — t/2, k! = 2mk, 7 = 271770, we get 

/ — 1 + e k'(x+py+u,t)+'y' + a 2 e -k'(x+py+ojt)-% + ^^'(z+py+wt)-^ 

Then the periodic solution (3.9) convergence to well-known soliton solution 

u = 2(\nf) xx , f = l + e fc '(^+-')+7o. (3.10) 
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Figure 9. One-periodic wave for KP equation and the effect of parameters on wave shape: 
(a) along x-axis, (b) along y-axis, (b) along t-axis, where k — 0.1, t = i, p = 2. 



We only need to prove that 

uj — > —4k 3 ir 2 + 3kp 2 , as a 



(3-11) 



In fact, it is easy to see that 

an = -32kn 2 (a 2 + Aa 8 + 
a 12 = 1 + 2a 2 + ■ ■ ■ , 

a 21 = -87r 2 £;(a + 9a 5 H ), 

a 22 = 2(a + a 5 + •••), 

h = 512A; 4 7r 4 (a 2 + 16a 8 + •••)- 9Qk 2 7r 2 p 2 (a 2 + 4a 8 + ■■■), 



b 2 = 32kV(a + 81a 5 + •••)- 24k 2 ir 2 p 2 (a + 9a b + ■■■). 
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Figure 10. One-periodic wave for KP equation and the effect of parameters on wave shape: 
(a) along x-axis, (b) along y-axis, (b) along t-axis, where k = 0.1, r = 3i, p = 1. 



Therefore we have 



^i«22 — b 2 ai2 = —32k tt a + 24tt k p a + o(a), an0 2 2 — 012021 = 8kir a + o(a), 



which lead to (3.11). 
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Figure 11. One-periodic wave for KP equation and the effect of parameters on wave shape: 
(a) along x-axis, (b) along y-axis, (c) along t-axis, where k = 0.2, r = i,p = 1. 



Now we consider two-periodic wave solution of the KP equation (N = 2). Let 
(fciPi, & = (fcik'i, ^2^2)- Substituting (3.4) into (3.3), we have 



Qj.j ^ ^ G(D D ^g27ri<^.n>+7ri<rn,n> _ g27ri<£,m>+7ri<Tm,m> 

m,n£Z 2 

= 2 G(2tt« < n - m, k >, 2m < n - m, p >, 2m < n - m, w >) e 27 ™<<>+ m >+ 7 "(<" Tl < m >+<™<"' 

m.n£Z 2 

00 

n+rn= m ' ^ ^ G(2ni < 2n - m' , k > ,2ni < 2n - m' , p > ,2ni < 2n - m' ,u> >) 

m'eZ 2 ni,n 2 =— 00 

x cxp(7r«(< r(n — m'), n — m' > + < rn, n >)) exp(27ri < £, m' >) 
= J] G(m / 1 ,n4)e 2,r<<£ ' m ' > =0. 
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It is easy to calculate that 

G(m' x ,m' 2 ) 

OO 

= ^ G ( 27Tl < 2fl ~ m '' k >' 2m < 2n - m '' P>,2™<2n- TO', W >) e Ti«r(n-m').n-"»'>+<™,n» 



111 ,712 — — oo 



= J2 G \ ^iJ2( 2n 'j - ( m 'j - 2<y,-i))fcj,27r*5^(2n$ - (m$ - 25 jl ))k jPj ,2m^{2n' :j - (m'j - 2S jl ))k J c 

ni,n2= — OO \ j — 1 j — 1 j — 1 

x exp < Tri [Oj- + 5ji)T jk {n' k + 4;) + ((m^ - - r^) + Sji)T jk {m' k - 2<5 fe( - n' fe ) + 4;)] 
[ i,k=i 

'G(mi - 2,m(,)e 27ri ( m 'i- 1 ) Tll+27rim 2 T i2, Z = 1 
^G(m[,m' 2 - 2) e 2 "( m 2-i)^2+2™m' 1 T 12; ; = 2 

which implies that if 

G(0, 0) = (5(0, 1) = (7(1, 0) = (5(1, 1) = 0, 
then G(m' 1 ,m' 2 ) = and (3.2) is an exact solution of the KP equation. Denote 

OO 

dji = —A-n 2 ki < 2n — m J , A; > (2ni — m{)5j(n), 

ni,n2=— oo 

OO 

Oj2 = — 47r 2 A; 2 < 2n — m J , /c > (2n 2 — m^t^n) 

ni,n2=— oo 
oo 

a j3 = 247r 2 ^ < 2n - m j , k > 2 Sj(n), 

ni,7i2=— oo 

OO 

%4 = ^ <$j(n) 
ni,ri2=— oo 
oo 

bj= {IQtt 4 <2n-m j ,k> 4 -12n 2 <2n-m j ,p> 2 )5 j {n), 



m,7i2=— oo 



^.(n) = e ^<-(«-™ 3 )" 3 >+^<™.™> 5 j = 1,2,3,4 
7^ = (0,0), m 2 = (l,0), m 3 = (0,l), m 4 = (1, 1) 
A = (a kj ) 4 x4, b = {biMMM) T , 
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then we have 



A 



L0 2 
W 



b, (3.12) 



from which we obtain 



Ax A 2 A 3 A 



= UJ2 = u = c = (3.13) 

where A = \A\ and Aj are produced from A by replacing its jth column with b 
respectively. Therefore, we arrive at two-periodic wave solution 

u = u + 2(\nf) xx , (3.14) 

where / and uj\,uj 2 are given by (3.4) and (3.13) respectively. The properties of the 
solution are shown in fig.12-fig.16. Fig. 12 and fig. 13 shown effect of parameters k\, k 2 
on period and shape of wave. Fig. 14 gives effect of parameter rn, ri 2 , t 22 on period and 
shape of wave. The fig. 15 and fig. 16 show that two-periodic wave can degenerate to 
one-periodic wave when k\ is sufficient small. 

Two-soliton solution of the KP equation can be obtained as limit of the periodic 
solution (3.14). We write / as 

j — l _|_ ( e 27 ™?l _|_ e -27ri§i^ e 7ririi _|_ ^2^2 _|_ e -2m^ e ^iT22 
_|_(g27ri(£i+6) _|_ e -27ri(€i+5 2 )^ e 7ri(r 11 +2ri2+r22) _)_... 

Setting 7^ = 7i + |rn, i 2 = 7 2 + |r 22 , ^ = 2m^i-mT n , £ 2 = 2iri£ 2 -7"r 22 , t 12 = it' 
(t' is a real), we get 

/ = 1 + e « + e £ + e £+&+2**l2 + a 2 e - 5l + a 2 e -^ + a 2 a 2 e - €l -6+2^n + . . . 

— ► 1 + + e«2 + eft-Ki-^', as ai ,a 2 — > 

where 

on = e" iT11 , a 2 = e niT22 , ^ = 2?ri (x + ^ + Ujt) + 7 •] , 
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Figure 12. Two-periodic wave for KP equation and the effect of parameters on wave shape: 
(a) along x-axis, (b) along y-axis, (c) along t-axis, where k\ = 0.2, ki — —0.3, 
Tn = 0.1i,Ti2 = 0.2i,r 22 = 3i,pi = l,p 2 = 2. 
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Figure 13. Two-periodic wave for KP equation and the effect of parameters on wave shape: 
(a) along x-axis, (b) along y-axis, (c) along t-axis, where k\ = 0.1, &2 = —0.3, 
Tn = 0.1i,Ti2 = 0.2i,r 22 = Zi,pi = l,p 2 = 2. 



References 

[1] M.J. Ablowitz and P.A. Clarkson, Solitons, Nonlinear Evolution Equations and Inverse Scattering, 

Cambridge University Press, 1991. 
[2] R. Beals and R.R. Coifman, Comm. pure. Appl. Math. 37 (1984) 39. 

[3] V.B. Matveev and M.A. Salle, Darboux transformation and solitons, Springer, Berlin, 1991. 

[4] C.H. Gu, H.S. Hu and Z.X. Zhou, Darboux Transformations in Soliton Theory and its Geometric 

Applications, Shanghai Science Technology Publisher, Shanghai, 1999. 
[5] S.B. Leble and N.V. Ustinov, J. Phys. A 26 (1993) 5007. 
[6] P.G. Esteevez, J. Math. Phys. 40 (1999) 1406. 
[7] R. Hirota and J. Satsuma, Prog. Thcor. Phys. 57 (1977) 797. 
[8] R. Hirota, Direct methods in soliton theory, Springer, 2004. 



22 



(a) 



(b) 




(c) 
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Figure 15. Two-periodic wave for KP equation and the effect of parameters on wave shape: 
(a) along x-axis, (b) along y-axis, (c) along t-axis, where k\ = 0.01, k 2 = —0.3, 
Tn = O.H,t 12 = 0.2z,r 22 =i,Pi = 0.2, p 2 = 0.3. 
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Figure 16. Two-periodic wave for KP equation and the effect of parameters on wave shape: 
(a) along x-axis, (b) along y-axis, (c) along t-axis, where k\ = 0.001, k 2 = —0.3, 
Tn = O.li, t\2 = 0.2i, r 22 =i,p\ = 0.2, p 2 = 0.3. 



